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I.  INTRODUCTION 


Flame  chemistry  models  for  even  simple  chemical  oxidations  include 
hundreds  of  chemical  reactions Commonly ,  rate  coefficient  data  is  available 
for  few  of  the  proposed  reactions.  In  some  cases,  x®_acti°n  network  schemes  do 
not  provide  even  qualitative  descriptions  of  the  observe3~chemist-ry-^-^,The 
methoxy  radical,  CH3O,  appears  as  a  constituent  in  chemical  models  that 
describe  the  oxidation  of  fuels,  especially  combustion  systems.  Numerous 
experimental^-^  and  theoretical^-^  research  projects  have  f-ocused  upon  this 


1.  K.  Iwasaki  and  K.  Tariyana,  "Electron  Spin  Resonance  of  Methoxy  Radicals 

Trapped  in  the  Methanal  Irradiated  at  4.2K,"  J ,  Am.  Chem.  Soc , ,  Vol.  100, 
p.  1964,  1978.  “ 

2.  D.K.  Russell  and  H. F.  Radford,  "Analysis  of  the  LMR  Spectra  of  Methoxy, 
CH20,"  J.  Chem.  Phys.,  Vol.  72,  p.  7250,  1980. 

3.  H.R.  Wendt  and  H.E.  Hunzicker,  "Electronic  Absorption  Spectrum  of  CH^O," 
J.  Chem.  Phys.,  Vol.  71,  p.  5202,  1979. 

4.  K.  Ohbayashi,  H.  Akimata,  and  I.  Tanaka,  "Emission  Spectra  of  CH^O,  C2H5O 
and  i-C^HyO  Radicals,"  J.  Phys.  Chem.,  Vol.  81,  p.  798,  1977. 

5.  M.  Sutoh,  N.  Washida,  H.  Akimoto,  M.  Nakamura,  and  M.  Okuda,  “The 
Emission  of  CH3O  from  Reactions  of  Metastable  Rare  Gas  Atoms  with 
Methanol,"  J,  Chem.  phys.,  Vol.  73,  p.  591,  1980. 

6.  G.  Imoue,  H.  Akimato,  and  M.  Okuda,  "Spectroscopy  of  the  CH3O  A^a^-X^e 
System  by  Laser-Excited  Fluorescence  Method,"  J.  Chem.  Ph/s.,  Vol.  72, 
p.  2769,  1980. 

7.  N.  Sanders,  J.E.  Butler,  L.R.  Pasternak,  and  J.R.  McDonald,  "CH3O 
Production  from  266  mm  Photolysis  of  Methyl  Nitrite  and  Reaction  with 
NO,"  Chem.  Phys.,  Vol.  49,  p.  17,  1980. 

8.  D.E.  Powers,  J.B.  Hopkins,  and  R.E.  Smalley,  "Laser  Production  of  Set- 
Cooled  Radicals.  Methoxy  and  Methoxy- Argon ,"  J.  Phys.  Chem.,  Vol.  85, 
p.  2711,  1981. 

9.  B.K.  Janousek,  A.H.  Zimmerman,  K.J.  Reed,  and  J. I.  Brauman,  "Electron 
Photodetachment  from  Aliphatic  Molecular  Anions.  Gas  Phase  Electron 
Affinities  of  Methoxyl,  Test-Butoxyl ,  and  Neopentoxyl  Radicals,"  J.  Am. 
Chem.  Soc.,  Vol.  100,  p.  6142,  1978. 

10.  P.C.  Engelking,  B.G.  Ellison,  and  W.C.  Lineberger,  "Laser  Photodetachment 
Electron  Spectrometry  of  Methoxide,  Deuteromethoxide  and  Thiomethoxide : 
Electron  Affinities  and  Vibrational  Structure  of  CH^O,  CD^O.a  nd  CH^S," 

J.  Chem.  Phys.,  Vol.  69,  p.  1829,  1978. 

11.  D.R.  Yarkony,  H.F.  Schaefer,  III,  and  S.  Rothenberg,  "Geometries  of  the 
Methoxy  Radical  (X^E  and  A^Aj  States)  and  the  Methoxide  Ton,"  J,  Am. 

Chem.  Soc. ,  Vol.  90,  p.  656,  1974. 
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12.  W.A.  Lathan,  C.A.  Curtiss,  W.J.  Hehre,  J.B.  Lisle,  and  J.A.  Pople, 
"Molecular  Orbital  Theory  of  the  Electronic  Structure  of  Organic 
Compounds,"  Progr.  Phys.  Org.  Chem. ,  Vol.  11,  p.  175,  1974. 

13.  H.  Uraeyaraa  and  S.  Nakagawa,  "Spectroscopic  Observation  of  the  Methoxy 
Radical,”  Chem,  Pharm,  Bull.,  Vol.  25,  p.  1671,  1977. 

14.  J.T.  Hougen,  "Double  Group  Considerations,  Jahn-Teller  Induced  Rovibroni 

Effects,  and  Nuclear  Spin-Electron  Spin  Hyperfine  Hamiltonian  for  a 
Molecule  of  Symmetry  C-iv  in  an  Electronic  ^E  State,"  J.  Mol.  Spectrosc. , 
Vol.  81,  p.  73,  1980.  y  ‘ 

15.  G.F.  Adams,  G.D.  Bent,  G.D.  Purvis,  and  R.J.  Bartlett,  "Calculation  of 
Dissociation  Energies  Using  Many- Body  Perturbation  Theory,"  Chem.  Phys. 
Letters ,  Vol.  81,  p.  461,  1981. 

16.  C.F.  Jaeckels,  "Electronic  Structure  Calculations  for  Methoxy:  The  X^E 
and  A^Aj  States,"  J.  Chem.  Phys.,  Vol.  76,  p.  505,  1982. 

17.  G.F.  Adams,  R.J.  Bartlett,  and  G.D.  Purvis,  "On  the  Unimolecular 
Reactions  of  CH^O  and  C^OH,"  Chem.  Phys.  Letters,  Vol.  87,  p.  311,  1982 

18.  G.D.  Bent,  G.F.  Adams,  R.H.  Bartram,  G.D.  Purvis,  and  R.J.  Bartlett, 
"Many-Body  Perturbation  Theory  Electronic  Structure  Calculations  for  the 
Methoxy  Radical.  I.  Determination  of  Jahn-Teller  Energy  Surfaces,  Spin 
Orbital  Splitting  and  Zeeman  Effect,"  J.  Chem.  Phys.,  Vol.  76,  p.  4144, 
1982. 

19.  S.  Saebo,  L.  Radom,  and  H.F.  Schaefer,  III,  "On  the  Weakly  Exothermic 
Rearrangement  of  Methoxy  Radical  (CH^OH),”  J.  Chem.  Phys.,  Vol.  78, 
p.  845,  1983. 
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radical.  Most  of  these  researches  attempted  to  elucidate  the  molecular 
structure  and  spectroscopy  of  this  radical;  tasks  which  remain  incomplete, 
however.  To  date,  no  empirically  derived  structure  has  been  obtained,  nor  has 
the  complete  vibrational  spectrum  been  determined.  While  such 
characterization  studies  continue,  the  chemical  reactions  that  include  the 
radical  populate  nearly  all  combustion  models.  Here  too,  the  lack  of 
experimental  and  theoretical  information  leads  to  great  uncertainty  in 
estimating  rate  coefficients  for  these  reactions.  Our  previous  research  on 
this  radical^2  was  motivated,  in  part,  by  uncertainty  in  rate  coefficient  data 
for  the  hydrogen  dissociation  reaction, 

ch3o  +  ch2o  +  H  (  1 ) 

Since  methoxy  is  but  one  hydrogen  atom  more  complex  than  formaldehyde,  we  also 
considered  the  1 ,2-hvdrogen  shift  reaction  that  isomerizes  methoxy  to  form 
hyd  roxymethy 1 

CH3O  ~  CH2OH  (2) 

As  was  the  case  with  the  formaldehyde  results, 2(2  the  electronic  energy  barrier 
for  methoxy  isomerization  wa  only  slightly  higher  than  the  dissociation 
barrier.  This  result  supported  Radford's  interpretation  of  experimental  data 
suggesting  that  the  isomerization  reaction  ( Eq .  2)  may  compete  with  the 
dissociation  reaction  ( Eq .  I).21  Subsequently,  Batt,  Burrows,  and  Robinson22 
estimated  an  activation  energy  for  the  isomerization  reaction,  and  concluded 
that  the  dissociation  reaction  dominates  at  high  temperature.  They  noted, 
however,  that  the  rearrangement  reaction  may  compete  with  the  dissociation  at 
high  temperature. 

This  report  describes  electronic  structure  calculations  pertinent  to 
these  two  reactions  of  methoxy.  The  quantum  chemical  results  are  used  to 
obtain  rate  coefficient  estimates  for  these  two  reactions  for  both  low  and 
high  temperatures,  in  the  low  pressure  limit  region. 


20.  .1 . D.  Goddard  and  H.F.  Schaefer,  III,  "The  Photodissociation  of 

Formaldehyde:  Potential  Energy  Surface  Features,"  J.  Chem.  Phys., 

Vo  1 .  70,  p.  51  17  ,  1979. 

21.  H.F.  Radford,  "The  Fast  Reaction  of  CH20H  with  02,”  Chem.  Phys.  Letters, 
Vol .  71,  p.  195,  1980. 

22.  L.  Batt,  J.P.  Burrows,  and  G.N.  Robinson,  "On  the  Isomerization  of  the 
Methoxy  Radical:  Relevance  to  Atmospheric  Chemistry  and  Combustion,” 
Chem.  Phys.  Letters,  Vol.  78,  p.  4b7 ,  1981. 


-v- 
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II.  THEORETICAL  CHEMISTRY  OVERVIEW 


The  description  of  chemical  reaction  pathways  using  quantum  chemical 
methods  has  been  simplified  by  the  introduction  of  analytic  formulas  for  the 
gradient  of  the  energy  and  the  second-derivatives  of  the  energy.  i  This 
research  employed  the  use  of  analytical  energy  gradients,  computed  with 
reference  to  unrestricted  Hartree-Fock  (UHF)  wavefunctions ,  to  locate  those 
points  on  the  energy  hypersurface  where  the  gradients  vanish.  Such  points 
correspond  to  equilibrium  structures,  when  there  are  no  negative  second- 
derivatives,  and  transition  state  structures,  when  there  is  one  and  only  one 
negative  second  derivative.  Structural  parameters  and  second-derivatives  for 
the  methoxy  radical,  and  the  two  transition  states  were  obtained  using  the 
GRADSCF  system  of  electronic  structure  codes. ^  These  UHF  calculations 
obtained  wavefunctions  expanded  in  a  6-3IG**  basis  set.~  Subsequent 
structural  optimizations  obtained  using  a  Dunning-Huzinaga  double-zeta  basis 
set^  augmented  by  polarization  functions  on  all  atoms^'  did  not  yield 
significantly  different  results. 

Self-consistent  field  methods,  like  UHF,  yield  unreliable  predictions  of 
chemical  energy  differences.  These  methods  exclude  the  effects  of  correlated 
electron  motions.  There  are  a  number  of  quantum  chemical  methods  available  to 
estimate  the  correlation-energy  contribution.  In  this  study,  we  have  used  the 
1 inked-cluster  based  method  known  as  many-body  perturbation  theory  (MBPT) . 

The  reference  function  for  the  MBPT  calculations  was  taken  to  be  a  UHF 
wavefunction  in  the  double-zeta  plus  polarization  basis  set.  The  MBPT 
calculations  are  fourth-order  calculations  that  include  all  single-,  double-, 
and  quadruple-excitation  diagrams  that  contribute  at  fourth-order. 


23.  M.J.  Page,  P.W.  Saxe,  B.H.  Lengsfield,  III,  and  G.F.  Adams,  "On  the 
Development  of  MCSCF  and  MCSCF-CI  Methods  for  the  Calculation  of  Hessian 
Matrices,"  to  be  published. 

24.  A.  Komornicki,  National  Resource  for  Computations  in  Chemistry  (NRCC) , 
Software  Catalog,  Program  Number  QH04  (GRADSCF),  Vol.  I,  1980. 

25.  P.C.  Hariharan  and  J.A.  Pople,  "Accuracy  of  AHj,,  Equilibrium  Geometries  by 
Single  Determinant  MO  Theory,"  Mol .  Phys. ,  Vol.  27,  p.  209,  1974. 

26.  T.H.  Dunning,  Jr.,  "Gaussian  Basis  Functions  for  Use  in  Molecular 
Calculations.  I.  Contraction  of  (9SSp)  Atomic  Basis  Sets  for  First-Row 
Atoms,"  J.  Chem.  Phys.,  Vol.  53,  p.  2823,  1970. 

27.  T.H.  Dunning  and  P.J.  Hay,  "Gaussian  Basis  Sets  for  Molecular 
Calculations  in  Modern  Theoretical  Chemistry,  Methods  of  Electronic 
Structure  Theory,"  Vol.  3,  H.F.  Schaeffer,  editor,  1977. 
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To  estimate  the  unimolecular  rate  coefficients  we  use  a  simplified  theory 
of  unimolecular  reactions  developed  by  Troe. 28  29  xroe's  review  paper^^  cites 
the  lack  of  detailed  knowledge  of  the  potential  energy  hypersurfaces  and  the 
intermolecular  energy  transfer  processes  as  deficiencies  that  hamper  the 
development  of  predictive  rate  coefficient  models.  The  electronic  structure 
calculations  reported  here,  and  in  references  15  and  17,  provide  useful 
information  about  the  potential  energy  hypersurfaces  for  the  methoxy 
radical.  Details  of  the  energy  transfer  processes  stand  beyond  the  scope  of 
this  work.  Troe's  model  avoids  the  energy  transfer  processes  by  invoking  the 
strong  collision  hypothesis,  and  seeking  ways  to  correct  this  overestimate. 


Troe^  expresses  the  limiting  low-pressure  pseudo-first-order  rate 
coefficient,  kg,  as 


k0  =  ec  C 


(3) 


sc 


where  $c  represents  the  temperature  dependent  collision  efficiency  and  k!! 
represents  the  strong-collision  form  of  the  limiting  rate  coefficient.  The 
strong-collision  rate  coefficient  is  expressed  as 


p(E  JkT 


koC  =  i*  zlj  ^ 


F_F  ,F  K  eXP  ("E  /kT>  > 
E  rot  anh  o 


(A) 


where  Z^j  is  the  Lennard-Jones  collision  frequency,  p(E0)  is  the  harmonic 
density  of  states,  Qv  is  the  vibrational  partition  function,  and  the  F  terms 
are  correction  factors  for  the  energy  dependence  of  the  density  of  states, 
rotational  effects  and  anharmonicity ,  respectively.  The  exponential  contains 
the  critical  energy,  EQ ,  the  energy  difference  between  the  reactant  and  the 
transition  state  measured  between  the  vibrational  zero-points.  In  this  work, 
we  used  expressions  given  by  Troe  in  Reference  27  to  evaluate  the  factors 
of  k;c. 

A  more  precise  approach  to  the  treatment  of  competing  reactions  requires 
that  we  treat  the  two  reactions  simultaneously.  Formally,  for  a  two-channel 
competing  reaction  process,  the  total  and  single-channel  rate  coefficients  are 
given  by^*-* 


28.  J.  Troe,  "Theory  of  Thermal  Unimolecular  Reactions  at  Low  Pressures. 
II.  Strong  Collision  Rate  Constants  Applications,"  J.  Phys,  Chem. , 
Vol .  66,  p.  4758,  1977.  . . 


29.  J.  Troe,  "Predictive  Possibilities  in  Unimolecular  Rate  Theory,"  J.  Phys. 
Chem. ,  Vol.  83,  p.  114,  1979.  ”  " 


30.  Th.  Just  and  J.  Troe,  "Theory  of  Two-Channel  Thermal  Unimolecular 

Reactions.  I.  General  Formulation , "  J.  Phys,  Chem.,  Vol.  87,  p.  3068, 
1980. 


(5) 


I 

!  ko  =  E/  ”  h( E)  f(E)  [kL(E)  +  k2(E)l  dE 


kio  ~  “  k20  ^ 

k20  =  p/ °°  h<E)  f(E)  k2(E)  dE  (7) 

02 

where  kjQ  represents  the  rate  coefficient  for  the  channel  with  the  lower 
critical  energy  and  k^Q  that  for  the  channel  with  the  higher  critical  energy, 
h( E)  is  the  population  of  reactive  states  above  the  lower  threshold  energy, 
and  f(E)  is  the  equilibrium  population  at  energy  E.  Specific  expressions  for 
the  total  rate-coefficient  are 


P(E  )kT 

ko  -  z  «*p  <-Eoi/kT>  fe  (dhr>  (8) 

where  Eq^  is  the  critical  energy  for  the  lower  critical  energy  channel, 
and  a  represents  the  average  energy  removed  in  down  transitions.  The 
corresponding  expression  for  the  upper-channel  rate  coefficient  is 


02 


=  exp  (- 


(Eo;-Eoi)<MFEkT) 
aF  kT 

C* 


■)* 


crt-F  kT 
[  ( - E — ^ 

1  aF  kT 
E 


k2(E) 

E02  <kl<E)*k2(E))* 


exp 


(E-Eno)(  «fF  kT) 

( - irr-L-)dEI 

h 


(9) 


The  integral  appearing  in  the  last  equation  can  be  approximated  in  various 
ways.  If  the  specific  branching  ratio 


k2(E) 

V2(E)  =  k  j"(  E)+k2"("E)'  » 


(10) 


varies  only  slightly  with  energy,  then  the  ratio  will  be  nearly  equal  to 
V2^E02^’  Vo(E)  varies  markedly  with  E,  then  the  specific  rate  constants 

kj(E)  and  k2(E)  must  be  used  in  evaluating  the  integral,  or  an  exponential 
model  must  be  used  to  approximate  i^CE).  Since  we’re  interested  in  the 
variation  of  the  rate  coefficients  with  temperature  over  a  large  temperature 
range,  the  evaluation  of  the  specific  branching  ratio  will  be  treated  in  a 
consistent  manner.  A  more  thorough  discussion  of  the  theory  of  two-channel 
unimolecular  reactions  is  given  by  Just  and  Troe.^  The  treatment  used  in 
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this  study  assumes  a  vibrational  two-channel  model  for  the  low-pressure  limit 
of  the  unimolecular  reaction.  Both  unimolecular  reactions  of  the  methoxy 
radical  are  characterized  by  rigid  transition  state  complexes. 


III.  ELECTRONIC  STRUCTURE  THEORY  RESULTS 

In  our  previous  report  on  the  unimolecular  reactions  of  the  methoxy 
radical,^'  we  reported  structural  parameters  for  CH^O  and  a  series  of 
transition  states  relevant  to  the  chemis  ry  of  the  radical.  As  noted  by 
Saebo,  et  al,^  the  structure  reported  in  Reference  17  for  the  methoxy  radical 
corresponds  to  the  ^A"  component  of  the  ^E  state,  rather  than  the  ^A'  ground 
state.  The  correct  structural  parameters  and  electronic  energy  are  reported 
for  the  X^A'  state  of  methoxy  in  this  report.  The  transition  state  structures 
cited  in  Reference  17  correspond  to  points  on  the  ^ A'  hypersurf ace .  Thus,  the 
electronic  contribution  to  the  critical  energies  equals  the  electronic  energy 
differences  between  the  two  transition  state  structures  and  the  equilibrium 
structure.  Table  1  contains  a  summary  of  the  structural  information  for  the 
radical  and  the  transition  states,  and  Table  2  contains  the  results  of  the 
electronic  energy  calculations  obtained  with  the  double  zeta  plus  polarization 
basis  set  at  several  levels  of  theory.  The  most  accurate  energy  differences 
should  be  those  predicted  by  fourth-order  MBPT;  thus,  the  electronic  energy 
barriers  equal  35.2  kcal/raole  and  35.6  kcal/mole  for  dissociation  and 
Isomerization,  respectively. 


TABLE  1.  OPTIMIZED  GEOMETRICAL  PARAMETERS  FOR  STATIONARY  POINTS 

ON  THE  METHOXY  SURFACE3 »b 


ch3o  (X2A’> 


rCC-Hj) 
r(C-H3) 
r( C-0) 
OCH1 
OCH3 

h1ch2 


1 .086 
1.089 
1.382 
111.7 
106.1 
110.5 


Transition  Structure:  Dissociation 


r(C-H1) 

1.087 

r(C-H3) 

1.834 

r( C-0) 

1 .225 

OCH1 

121.4 

OCH3 

99.8 

Hj  CH2 

116.7 

Transition  Structure :  Isomerization 


rCC-Hj) 

r(C-H3) 

r(C-0> 

0CH1 

OCH3 

Hj  CH2 


1.078 

1.264 

1.367 

117.4 

53.5 

118.3 


-  Units: 
a  -  Units: 


Bond  Lengths,  Angstroms. 
Bond  Angles,  Degrees. 
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TABLE  2.  ELECTRONIC  ENERGIES  (UNITS:  HARTREES) 


Transition 

Transition 

Radical 

State 

State 

Equilibrium 

Dissociation 

Isomerization 

UHF 

-114.4484 

-114.3845 

-114.3668 

SDQ-MBPT(4) 

-114.7621 

-114.7060 

-114.7054 

To  determine  the  critical  energy  for  reaction  requires  knowledge  of  the 
zero-point  vibrational  energy  for  the  reactant  and  transition  state.  By 
computing  the  second-derivatives  of  the  molecular  energy  with  respect  to 
nuclear  distortion,  we  can  determine  the  harmonic  vibrational  frequencies  for 
these  structures.  The  harmonic  frequencies  for  the  methoxy  radical  and  the 
two  transition  states  are  presented  in  Table  3.  Note  that  each  of  the 
transition  states  has  one,  and  only  one,  imaginary  frequency.  With  these 
results,  we  estimate  the  critical  energy  for  dissociation  to  be  29.5 
kcal/mole,  while  that  for  isomerization  is  31.7  kcal/mole. 


TABLE  3.  VIBRATIONAL  FREQUENCIES  FOR  POINTS  ON  THE  METHOXY  HYPERSURFACE 

(UNITS:  WAVENUMBERS) 


Equilibrium 


Dissociation  TS* 


Isomerization  TS* 


3274 

3768 

3366 

3255 

3138 

3249 

3188 

1708 

2473 

1668 

1518 

1621 

1604 

1345 

1266 

1585 

1211 

1241 

1283 

668 

1074 

1226 

504 

945 

1130 

955i 

25691 

transition  State. 


To  anticipate  our  treatment  of  the  rate  coefficients,  we  note  that  the 
vibrational  frequencies  for  the  methoxy  radical  share  the  defect  of  all 
frequencies  determined  using  SCF  analytical  gradients.  Such  frequencies  are 
estimates  of  harmonic  vibrations,  while  empirically  observed  vibrational 
frequencies  include  anharmonic  character.  The  experience  of  quantum  chemists 
in  computing  frequencies  suggests  that  the  results  reported  in  Table  3  will 
exceed  observed  frequencies  by  approximately  10%.  Therefore,  we  reduce  the 
predicted  frequencies  before  use  in  the  rate  coefficient  estimation. 


IV.  RATE  COEFFICIENT  ESTIMATES 


Sy»: 

,*-\v 

• ' 

PT"1 


The  rate  coefficients  for  dissociation  and  isomerization  of  the  methoxy 
radical  for  temperatures  in  the  range,  300K  <  T  <;  1200K,  are  evaluated  in  two 
ways.  We  present  first  the  results  of  separate  calculations  of  the  two 
processes,  ignoring  the  complications  introduced  by  formal  consideration  of 
the  competitive  nature  of  the  reactions.  We  then  consider  the  reactions  as 
two  competitive  reactions,  determining  branching  ratios  using  Troe's 
vibrational  two-channel  model.  Comparison  of  the  two  treatments  allows  us  to 
assess  the  utility  of  the  simple  mode!  as  a  quick  and  dirty  tool  for  assessing 
branching  ratios. 

In  addition,  we  consider  also  the  effect  upon  these  results  of  errors  in 
the  quantum  chemical  calculations.  It  is  commonlv  accepted  that  electronic 
structure  methods  overestimate  barriers  to  unimolecular  reactions.  Thus,  it 
is  necessary  to  determine  whether  the  reasonable  correction  of  the  quantum 
chemical  barrier  heights  leads  to  significantly  different  relative  reaction 
rates,  or  branching  ratios. 

There  occurs  in  the  literature  several  values  for  the  heat  of  formation  of 
the  methoxy  radical.^  ^1  These  values,  when  combined  with  the  heats  of 
formation  of  formaldehyde  and  hydrogen  provide  the  heat  of  reaction  for  the 
dissociation.  Using  the  various  experimental  values,  the  heat  of  reaction  may 
be  22.2  kcal/mole  or  26.7  kcal/mole.  In  a  paper  on  the  a  priori  determination 
of  dissociation  energies,  we  estimated  the  heat  of  formation  of  CH^O  to  be 
1.0  ±  3.0  kcal/mole.^  This  value  has  the  virtue  that  the  error  bars 
encompass  the  experimental  values;  this  leads  to  a  dissociation  energy 
estimate  equal  to  25.1  kcal/mole.  Correcting  this  value  for  classical 
translational  energy  contributions,  we  obtain  the  OK  dissociation  energy, 

24.2  kcal/mole.  This  sets  a  lower  limit  on  the  critical  energy  for  this 
reaction  channel.  The  value  predicted  by  the  electronic  structure 
calculations  for  the  critical  energy  corresponds  to  a  recombination  barrier  of 
5.3  kcal/mole,  a  value  several  kcal  higher  than  the  barrier  estimated  for  the 
isoelectronic  ethyl  radical  barrier.  Thus,  there  is  some  evidence  that  the 
critical  energy  for  dissociation  has  been  overestimated  by  several  kcal/mole. 

Estimating  the  error  for  the  isomerization  reaction  barrier  is  less 
straightforward.  Although  we  expect  that  the  value  exceeds  the  "true" 
critical  energy,  it  is  difficult  to  make  a  reasonable  estimate  of  the  error. 

On  the  one  hand,  we  believe  that  the  theoretical  methods  used  provide  a  more 
reliable  prediction  of  the  isomerization  barrier  than  for  the  dissociation 
barrier.  Conversely,  we  have  not  included  the  contribution  of  triple 
excitations  on  the  fourth-order  energy  contribution;  an  effect  that  should 
reduce  the  isomerization  barrier  more  than  it  reduces  the  dissociation 
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barrier.  In  any  case,  however,  we  expect  both  barriers  to  be  reduced  by 
better  theoretical  calculations.  Therefore,  to  provide  a  useful  test  of  the 
effect  of  quantum  chemical  errors  on  the  branching  ratios,  we  assume  that  the 
''correct’’  critical  energies  are:  dissociation  27.0  kcal/mole;  isomerization, 
30.0  kcal/mole.  These  values  will  be  used  to  estimate  branching  ratios  for 
T  =  I200K. 

In  the  following  sections,  we  present  the  results  obtained  for  the  rate 
coefficients  by  treating  the  reactions  as  independent  and  coupled. 

A .  Independent  Rate  Coefficient  Treatment 

H  C 

The  data  used  to  evaluate  the  factors  in  the  expression  for  k  were 
determined  entirely  from  the  results  of  the  quantum  chemical  calculations 
which  are  summarized  In  Tables  1-3.  To  estimate  a  value  for  the  collisions 
efficiency,  we  assumed  that  the  average  energy  transfer  in  a  down  collision 
was  4.6  kj/mole  (1.099  kcal/mole),  and  used  Troe's  relationship  connecting 
&c  and  a.  We  summarize  in  Tables  4  and  5  the  factors  obtained  for  the 
dissociation  and  isomerization  reactions,  respectively.  The  values  of  k^, 
given  in  the  right-hand  column,  lead  to  the  conclusion  that  the  dissociation 
dominates  at  room  temperature,  but  that  the  isomerization  reaction  competes 
with  dissociation  at  high  temperatures.  To  some  extent  the  result  parallels 
the  change  in  relative  values  of  the  exponential  factor,  but  one  should  note 
that  the  harmonic  density  of  states  factor  always  favors  the  isomerization 
reaction.  In  fact,  the  density  of  states  factor  will  aLways  favor  the  channel 
with  the  higher  critical  energy,  since  we've  evaluated  the  density  of  states 
at  E0.  To  recapitulate,  this  model  predicts  that  the  isomerization  and 
dissociation  reactions  compete  at  high  temperature,  a  result  that  concurs  with 
Batt’s  conjecture. 

What  is  the  effect  on  the  high-temperature  rate  coefficients  of  reducing 
the  critical  energies  for  the  two  reactions?  The  detailed  results  for  the  two 
rear.;  ms  are  collected  in  Table  6.  The  ratios  of  the  k.()’s  changes  little 
from  that  obtained  using  the  original  quantum  chemical  critical  energies. 
Apparently  the  results  of  assuming  a  greater  separation  between  the  two 
critical  energies  cancel  the  results  of  reducing  these  values.  Thus,  ii  we 
assume  that  the  branching  ratio  for  the  competitive  reactions  can  be  estimated 
by  computing  the  rate  coefficients  independently,  we  conclude  that  the 
isomerization  and  dissociation  reactions  compete  at  high  temperature. 

B.  Vibrational  Two-Channel  Model 


The  expression  for  the  total  unimolecular  rate  coefticient,  Kq .  (H),  is 
very  similar  to  the  expression  derived  by  Troe  for  the  low  pressure  rate 
coefficient  of  a  single  unimolecular  reaction.  Missing  are  rotational  and 
anharmonic  corrections.  The  expression,  however,  contains  a  term  relating  to 
the  ene rgy— t rans f er  step  size.  Since  we  estimated  this  value  to  determine  the 
collision  efficiency,  6C,  consistency  demands  retaining  the  same  value  in  this 
treatment . 

To  evaluate  the  expression7  for  the  rate  coefficient,  k™,  requires 
evaluating  the  specific  branching  ratio,  V2(E),  given  by  Eq.  1 10).  The 
expression  for  V2(E),  may  be  rewritten  as 


2V  J 


W1(E)+W2(E)  ’ 


vn/ 


where  Wj(E)  indicates  the  number  of  open  vibrational  reaction  channels  for  the 
products  of  the  lower  critical  energy  pathway.  Thus,  one  can  determine  V2(E) 
by  counting  the  number  of  open  channels  at  each  energy  for  each  pathway. 

Since  the  critical  energies  for  the  two  reaction  paths  are  so  nearly  equal, 
the  specific  branching  ratio  should  vary  only  slightly  with  energy  near  Eq2» 
so  we  follow  Troe  and  assume  that 


v2(e)  =  v2(e02) 

Finally,  rather  than  count  the  number  of  reactant  channels  open  for  each 
pathway,  we  assume  that  the  specific  branching  ratio  can  be  approximated 
substituting  the  vibrational  partition  functions  for  the  two  activated 
complexes.  This  approximation  has  no  formal  justification,  but  since  we  have 
estimates  of  the  vibrational  eigenvalues  orthogonal  to  the  reaction 
coordinate,  we  assume  that  the  populations  of  these  states  provides  an 
approximation  to  the  number  of  available  channels  for  each  pathway.  This 
leaves  us  to  evaluate  the  trivial  integral, 
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exp 
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EkT 


(12) 


Evaluating  the  expressions,  using  the  data  from  Table  5,  we  find  that  for 
T  =  300 


k2Q  a  .002  kQ  , 
while  for  T  =  1200K 


k2Q  —  .023  kQ 


Since 


k10  “  ^  “  k20  » 

it  is  obvious  that  the  two-channel  model  does  not  predict  the  isomerization 
reaction  to  compete  with  the  dissociation  reaction,  even  at  T  =  1200K.  This 
result  depends,  however,  upon  the  estimates  of  the  energy  transfer  parameter, 
a,  the  dilemma  faced  when  trying  to  use  unimolecular  models  as  predictive 
models . 

The  electronic  structure  calculations  summarized  in  this  report  provide 
information  about  the  potential  energy  hypersurface  for  the  methoxy  radical 
reactions  that  has  not  been  available  heretofore.  To  some  extent,  Troe's 
desire  for  details  of  the  potential  energy  surface  has  been  satisfied. 
Nonetheless,  two  simplified  models  for  rate  coefficient  estimates  provide 
qualitatively  different  descriptions  of  the  relative  importance  of  the  two 
reaction  channels.  This  result  demonstrates  the  need  for  more  p, ocise 
descriptions  of  the  energy  transfer  processes  that  dominate  the  low-pressure 
unimolecular  reactions. 


V.  CONCLUSIONS 


<  - 
/ 

The  results  of  this  study  do  not  pro  ide  a  definitive  demonstration  that 
the  isomerization  reaction  of  the  methoxy  radical  competes  with  the 
dissociation  reaction.  One  simplified  dynamical  model  suggests  that  the 
reactions  do  compete  at  high  temperature.  However,  Troe's  model  for  two- 
channel  thermal  unimolecular  reactions  predicts  that  the  isomerization 
reaction  does  not  compete  with  the  dissociation  reaction.  The  latter  result, 
however,  depends  strongly  upon  the  details  of  the  energy-transfer  processes 
that  control  the  low-pressure  rates,  and  on  the  rather  gross  approximations 
we've  used  in  applying  that  model. 

Variation  of  the  predicted  critical  energies  to  values  that  intuition 
suggests  are  closer  to  reality,  leads  to  virtually  no  change  in  the  branching 
ratio  predicted  by  the  uncoupled  dynamical  model.  The  rate  coefficient 
estimates,  however,  do  change  significantly,  indicating  that  predictive  rate 
coefficient  models  yet  await  more  accurate  quantum  chemical  calculations. 
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